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ABSTRACT 

The gas at the surfaces of molecular clouds in galaxies is heated and dissociated by photons 
from young stars both near and far. H i resulting from the dissociation of molecular hydrogen 
H2 emits hyperfine line emission at 21 cm, and warmed CO emits dipole rotational lines such as 
the 2.6 mm line of CO(l-O). We use previously developed models for photodissociation regions 
(PDRs) to compute the intensities of these H i and CO(l-O) lines as a function of the total volume 
density n in the cloud and the far ultraviolet flux Gq incident upon it and present the results in 
units familiar to observers. The intensities of these two lines behave differently with changing 
physical conditions in the PDR, and, taken together, the two lines can provide a ground-based 
radio astronomy diagnostic for determining n and Gq separately in distant molecular clouds. This 
diagnostic is particularly useful in the range Go < 100, 10 cm~^< n < 10^ cm"'^, which applies 
to a large fraction of the volume of the interstellar medium in galaxies. If the molecular cloud 
is located near discrete sources of far-UV (FUV) emission, the PDR-generated Hi and CO(l-O) 
emission on the cloud surface can be more easily identified, appearing as layered "blankets" or 
"blisters" on the side of the cloud nearest to the FUV source. As an illustration, we consider 
the Galactic object G216 -2.5, i.e. "Maddalena's Cloud", which has been previously identified 
as a large PDR in the Galaxy. We determine that this cloud has n « 200 cm^"^ and Gq « 0.8, 
consistent with other data. 

Subject headings: atomic processes - galaxies: ISM - ISM: atoms - ISM: clouds - ISM: molecules - 
molecular processes 
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1. Introduction 

The interstellar medium (ISM) in galaxies is 
excited, dissociated, and ionized by far-ultraviolet 
(FUV) photons produced by young O and B stars. 
Atomic gas in the ISM recombines into molecular 
form mainly through the catalytic action of dust 
grain surfaces. In particular, hydrogen nuclei in 
the ISM cycle repeatedly from the molecular (H2) 
to the atomic (H i) phase and back again, at rates 
depending on the incident FUV flux Go, the total 
volume density n of the gas, and the dust-to-gas 
ratio S. Regions in the ISM where the physics is 
dominated by FUV photons are called photodisso- 
ciation regions (PDRs). The surfaces of giant mo- 
lecular clouds (GMCs) are important (and ubiq- 
uitous) examples of PDRs in galaxies. 

The physics of PDRs has been explored in de- 
tail over the past ^ 30 years by many workers 
including D. J. HoUenbach, B. T. Draine, A. Dal- 
garno, J. H. Black, A. G. G. M. Tielens, E. van 
Dishoeck, J. Le Bourlot, and their students and 
collaborators. One major focus of this work has 
been to explain the ~ 1% line-to-continuum ra- 
tios of the far-infrared lines of [CII] and [Ol] in 
Galactic sources observed from high-altitude air- 
craft and balloon platforms. With the advent of 
more sensitive space-based observations, the mod- 
els were extended to include the weak rotational- 
vibrational spectrum of excited H2 observed on the 
active PDR surfaces of GMCs that are exposed to 
relatively intense FUV fluxes from nearby young 
stars. An excellent review of the observational and 
theoretical state of the field is given by Hollenbach 
& Tielens (1999). 

PDR model computations have become very de- 
tailed and comprehensive, raising the possibility 
that observations of multiple spectral lines includ- 
ing those from trace molecules such as CO can be 
used to "invert" the models in order to determine 
physical conditions in the ISM. For example, a set 
of such models has been computed by Kaufman 
et al. (1999) for use in the interpretation of far- 
infrared and submillimeter spectra from Galactic 
and extragalactic sources. However, we wish to 
point out that, although there is qualitative agree- 
ment on most of the physics, there are differences 
in the details of the various model codes currently 
in use around the world as well as differences in 
the numerical values of the parameters adopted. 



These differences can lead to disagreements in the 
values derived for physical conditions even for the 
same observational data: an example of such dif- 
ferences is illustrated in Li et al. (2002). Efforts 
are underway to intercompare models from differ- 
ent groups in a more systematic manner (E. van 
Dishoeck 2003, private communication). The re- 
sults obtained from any one model computation 
(such as the one we use here) can therefore be re- 
vealing in general, but the specific numerical val- 
ues obtained should be taken with some caution. 

The purpose of the present paper is to draw 
attention to the fact that the Hi produced from 
photodissociated H2 can also be a useful diagnos- 
tic for determining physical conditions in PDRs. 
In particular, the combination of the 21 cm Hi 
and 2.6 mm CO(l-O) lines provides a means of in- 
dependently estimating n and Go in distant PDRs 
using ground-based radio astronomy data. These 
two radio lines form a diagnostic that is particu- 
larly useful in the area of parameter space where 
Go < 100, 10 cm-3< n < 10^ cm-^; this range 
of FUV flux and total density is representative of 
most of the volume of the ISM in galaxies. 

2. PDR Models 

Our point of departure is the paper by Kaufman 
et al. (1999) that describes the results of steady 
state computations on a specific, extensive set of 
PDR models over a wide range of physical condi- 
tions. These models are simple one-dimensional 
semi-infinite slabs of gas with constant density 
n = n(Hi) + 2n(H2) of H nuclei subjected to an 
equivalent one dimensional flux of FUV photons 
Go measured in units of 1.6 x 10~^ ergs cm~^ 
over the photon energy range 6-13.6 eV.^ The line 
brightness is calculated along a ray normal to the 
surface of the slab. See Kaufman et al. (1999) for 
further details of their models. 

2.1. 2.6 mm CO(l-O) Line Emission 

The incident FUV flux dissociates and heats 
the cloud surface. At some depth into the cloud, 
the dissociation rate has decreased enough to al- 
low CO molecules to survive, and the density and 

^In these units, the average interstellaj: radiation field over 
47r sr is Go 1.7 (Draine (1978), see our Appendix B and 
Draine's Fig. 3). 
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heating rate arc still high enough to excite rota- 
tional line emission from these molecules. Among 
many important results, Kaufman et al. (1999) 
have computed the 2.6 mm C0(1 0) brightness, 
Ico, emanating from the cloud surface; this is the 
most widely studied molecular line in the ISM. 
In Figure la wc present their results for /co as 
a contour diagram of constant surface brightness 
(see their Fig. 11); we have extended the range of 
n and Gq down to 1 and 0.03, respectively, and 
labeled the contours with their original "theoret- 
ical" units of ergs cm~^ sr~^, as well as with 
"CO observer" units of K km s~^. This latter con- 
version is made according to: 

\3 

T^b X Ay = ;^/co = 6.4 x 10« x 7co, (1) 

2kb 

where Tmb is the measured emergent Rayleigh- 
Jeans brightness temperature assuming the cloud 
is completely resolved, AV is the profile velocity 
width in km s~^, A = 2.60 mm, and /cb is Boltz- 
mann's constant. 

We note that Figure 1 assumes that the dust- 
to-gas ratio S is that of the local ISM near the 
Sun, and only one side of the slab contributes to 
the observed emission. A further caution concerns 
the value of the density n inferred from a specific 
observed value of the CO(l-O) intensity in Figure 
la. There is a general trend in most PDR mod- 
els to underestimate the CO abundance, so that 
the density inferred from a given observed value 
of the CO (1-0) line is often higher than that ob- 
tained from other independent diagnostics. This 
is another example in which comparisons among 
the different modeling efforts will be very useful. 

The CO(l-O) surface brightness in Figure la 
is governed by a complex interplay of processes 
involving the dissociation and reformation of CO 
molecules, heating and cooling in the ISM, and 
the detailed excitation and macroscopic radiative 
transfer of the CO(l-O) hue. By contrast, the 
21 cm Hi line intensity is somewhat easier to com- 
pute, since in general this line is optically thin and 
independent of the kinetic temperature of the ISM; 
we can therefore ignore the physics of the 21 cm 
line formation and simply compute the amount of 
Hi produced. 



2.2. 21 cm H I Line Emission 

Incident FUV photons in the wavelength range 
912-1108 A (13.6-11.2 eV) also dissociate the H2 
on the surface of the GMC.^ This process was first 
described in the literature by Stecher & Williams 
(1967) and is by now well known: it has been sum- 
marized by Hollenbach & Tielens (1999), who also 
provide several other references. The dissociated 
H2 creates a "blanket" of atomic gas on the sur- 
face of the GMC; the Hi atoms recombine into 
H2 on dust grains and are then returned to the 
gaseous ISM. Approximate analytic expressions 
for the steady state column density N{R i) of the 
atomic gas have been derived by several authors 
under a variety of limiting assumptions (e.g. Stern- 
berg (1988); Goldschmidt & Sternberg (1995)) and 
are discussed in more detail in Appendix A. Com- 
putations of iV(Hl) in the context of the present 
"standard model" were first presented in Wolfire 
et al. (1990) (see their Fig. 6). An updated cal- 
culation using the latest code and the model pa- 
rameters in Kaufman et al. (1999) is presented as 
the dotted curves in Figure 3. The results differ 
slightly from those of Wolfire et al. (1990), owing 
to small differences in the assumed values of dust 
abundance, atomic abundances, and a different H2 
formation rate. This leads to changes in the depth 
of the Hi - H2 transition, and hence to different 
values of iV(Hi); the differences are, however, gen- 
erally less than about a factor of 2. 

A problem arises in attempting to combine the 
computational results for CO (1-0) with those for 
Hi. The "standard model" used here is bathed 
in a uniform flux of cosmic rays producing a con- 
stant ionization rate^ of 1.8 x 10"^'' s~^. These 
cosmic rays also produce a faint extended distri- 
bution of H I atoms owing mainly to ion chemical 
reactions everywhere in the cloud. However, as wc 
see below, such a faint extended distribution of H i 
will in general be removed from the observations 
by the data reduction process, so we would pre- 
fer to use models for A''(Hi) without a cosmic-ray 
component. Nevertheless, cosmic-ray ionization is 
important for the chemistry in the model, thereby 

^Photons with energies as low as 6.6 eV can also contribute 
to the dissociation at high FUV flux levels of Go > 10* 
(ShuU 1978). 

^This parameter is unfortunately missing from Table 1 of 
Kaufman et al. (1999). 
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controlling to some extent the abundance of CO 
molecules, so we ought to keep it in the computa- 
tions of the CO(l-O) line. 

As a compromise we have chosen to use the ap- 
proximate analytic expression for iV(Hl) since it 
does not include any effects of cosmic rays. This 
approximation is presented in Appendix A where 
it is compared with the standard computational 
model: it is an excellent fit to the model over the 
full range of parameter space important for this 
paper and provides a better context in which to 
compare the theory with the observations. The 
adopted result for the H i column is shown in Fig- 
ure lb. 

The conversion from H i column density in Fig- 
ure lb to 21 cm emission intensity in K km 
is: 

T,„b X Ay = 5.5 X 10"^^ X iV(Hi), (2) 

for A^(Hi) in units of atoms cm~^, and assuming 
the optical depth of the atomic gas to the 21 cm 
line radiation is small. 

A noteworthy feature of Figure lb is the con- 
stancy of N(H\) for a constant ratio of Gq/u (see 
also cq. A2); the same Hi column density can be 
produced in high FUV flux, high-density environ- 
ments or in low FUV flux, low-density environ- 
ments. In addition, the Hi column decreases with 
increasing density at a given value of FUV flux: 
this occurs because the destruction rate for H2 
varies as n, but the formation rate varies as 
(since the grain density is proportional to n). 

2.3. A Combined Model 

Figures la and lb show that the contours of 
constant CO(l-O) emission Iqq and constant Hi 
column density A'^(H l) behave in a complementary 
fashion. This suggests that over much of the n- 
Go plane an observational determination of Ico 
and A'^(Hl) for the same PDR would permit inde- 
pendent estimates to be made of both n and Go . 
In Figure 2 we show the two model computations 
from Figures la and b superposed, along with sev- 
eral additional contours and a "box," which will 
be discussed further below. 

2.3.1. Range of Validity of the Models 

The models become progressively less reliable in 
the top left corner of the plots for both Iqo and 
A''(Hi). This occurs because the assumption of 



a steady state model becomes increasingly unreli- 
able in the high Go - low n section of the diagram; 
in this region, radiation pressure causes grains to 
drift with speeds of the same order as the gas tur- 
bulence speeds (see Kaufman et al. (1999), Fig. 
2). Furthermore, this is also the area of the dia- 
gram where large H i column densities are predict- 
ed; however, values of iV(Hi) in excess of a few 
xlO^^ cm~^ are not likely to be observed since 
the 21 cm line may become optically thick in this 
regime: 

/oo 
TsT{v)dv 
-00 

1.82 X IO^^T^tAv 
« a few X 10^^ cm-^ 

for spin temperatures of Tg f« 50 — 100 K, profile 
FWHMs of Av « 4-8 kms-\ and r w 2-3, val- 
ues that probably represent the extremes typical 
for the observations. 

2.3.2. Resolution Effects 

It should be emphasized that the present calcu- 
lations of surface brightness can be applied direct- 
ly only to observations for which the sources of 
emission are resolved. In any realistic situation, 
the emission probably arises from several PDR 
surfaces. In this case, the density and FUV in- 
tensity derived from our model using observations 
of Ico and A'' (Hi) represent ensemble averages of 
the distribution of values for n and Gq. Wolfirc 
et al. (1990) have shown that these averages are 
somewhat biased toward higher FUV fields; bias 
toward high or low density depends on the cloud 
distribution. 

2.3.3. Interpretation of n 

The X axis of Figure 2 refers to the total volume 
density in our computational model of an isopyc- 
nic cloud. In a real GMC this is likely to be close 
to the actual density at the depth in the cloud 
where the CO(l-O) line is formed. However, it is 
not likely to be representative of the density of dis- 
sociated Hi. This is because the various heating 
processes at the surface of the GMC (including 
photodissociation) will lead to gradients in tem- 
perature and, therefore, if the pressure is roughly 
constant, to gradients in density. For example, the 
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model computations of Allen et al. (1995) for rela- 
tively low FUV fluxes^ of Go ~ 0.1 show that the 
density at the CO(l-0)-emitting layer is about 5 
times that at the H i - H2 interface. Their Figure 
4 shows that the CO emission emanates from a 
depth such that Ay > 2, where the temperature 
is just below 5 K (sec their Fig. 3b), whereas their 
Figure 2 shows the H i - H2 interface is effectively 
at very small values of Ay » 10~^, where the tem- 
perature is about 25 K. We have calculated tem- 
peratures for higher FUV flux levels Go = 3: the 
surface temperatures are higher in this case, » 60 
K at the Hi - H2 interface, and « 7.5 K at the 
CO emission layer, for a ratio of « 8. This ratio 
will continue to rise for yet higher photon fluxes, 
since the increased surface heating will result in 
further increases in surface temperature, but the 
increase in temperature deeper in the clouds where 
the C0(1 0) emission is formed will be less rapid 
because of more effective cooling. If the density 
and temperature are linked by an isobaric condi- 
tion, the density observed in the bulk of the Hi 
may be a factor of 5-10 lower than that inferred 
for the CO(l-O) line. 

It is clear from this discussion that constant- 
density models for PDRs such as the one we have 
assumed here are not likely to be very realistic. 
More sophisticated models, for example assuming 
pressure equilibrium, or a full computation with 
various forms of micro- and macro-turbulent pres- 
sure (e.g., Wolfire et al. (1993)) over the range of 
the Go-n parameter space may provide more re- 
alistic results and ought to be explored. 

3. Discussion 

3.1. Complementarity 

Figure 2 clearly shows the complementarity of 
21 cm H I and 2.6 mm CO(l-O) radio line emission 
as diagnostics of PDRs. Over the range of valid- 
ity of our calculations, the CO(l-O) line bright- 
ness Ico depends mostly on the volume density of 
the gas for n < 10^ cm~^, no matter how intense 

*The symbol x used by Allen et al. to denote the FUV flux 
has very nearly the same value as we use here, x = Go/0.85 
for irradiation over 27r sr. See Appendix B for a discussion 
of the various definitions of FUV values used in the litera- 
ture on PDRs. 



the FUV flux."'' On the other hand, for constant 
n the H i brightness increases nearly linearly with 
increasing FUV fluxes up to Go » 10—100, becom- 
ing logarithmic for higher FUV fluxes. According- 
ly, the combination of measurements of A^(H i) and 
7co from the same PDR can provide unique val- 
ues for n and Go over a substantial range in the 
diagram. 

3.2. Identifying the Relevant Hi and CO 

Emission 

How can extraneous emission from gas that lies 
far outside the cloud be identified? The spatially 
layered structure of PDRs provides a direct clue, 
as has been convincingly shown e.g. by the obser- 
vations of the Orion Bar region (see Fig. 2 in Hol- 
lenbach & Tielens (1999)) on the ^ 1 pc scale. A 
more or less edge-on viewing angle is required for a 
certain identification, and one ought to remove any 
extended emission that is associated with other 
clouds along the line of sight or produced by other 
excitation mechanisms. Galaxies that are viewed 
at low to intermediate inclination angles arc fa- 
vorably oriented for this identification to work, 
and Allen et al. (1997) and Smith et al. (2000) 
have shown that the signature morphology can be 
found at 100 pc scales in the highest resolution 
Hi images of nearby galaxies M81 and MlOl. As 
to the CO, it is generally assumed that this arises 
in PDRs, although the morphological signature is 
not always clear. 

3.3. G216 -2.5: An Example 

Figure 2 provides a method for determining n 
and Go for a specific PDR if iV(Hl) and Ico can 
be determined obscrvationally. As an example, we 
consider the Galactic object G216 -2.5 (also called 
"Maddalena's Cloud") which has been proposed 
to be a large PDR in the Galaxy by Williams & 
Maddalena (1996). 

3.3.1. Hi Columm, Density 

Determining the column density of that part of 
the Hi associated with G216 -2.5 is made quite 

^The inscnsitivity of Iqo to the FUV flux comes about be- 
cause increasing Go merely leads to more photodissociation 
of the surface CO, and the CO-emitting layer retreats deep- 
er into the static cloud where the temperatures are about 
the same as they were for the case of lower Go. 
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difficult by confusion from unassociatcd H I super- 
posed along the line of sight through the Galaxy. 
Several estimates can be made from the single-dish 
observations reported by Williams & Maddalena 
(1996). First, and perhaps simplest, we estimate 
that fraction of the total averaged H i profile shown 
in their Figure 1 that is in the velocity range of the 
CO (1-0) emission (also shown on the same fig- 
ure); the Hi profile shows three overlapping peaks 
and integrates to « 48 K x 58 km s-i« 2800 K 
km s~^= 5.1 X 10^^ cm~^. The averaged Ico Pro- 
file corresponds to the central Hi peak and is 8.5 
km s~^ wide (Maddalena & Thaddeus 1985), so 
the average amount of H i column associated with 
the CO emission is, by this estimate, probably not 
more than « (8.5/58) x 5.1 x lO^i = 7.5 x lO^o 
cm-2. In their §3.2.2, Wilhams & Maddalena 
(1996) carry out a more involved analysis with 
area integrations of the H i in the velocity range 
16-38 km s~^ and conclude that the "excess" Hi 
associated with the PDR is » 2 x 10^° cm"^; the 
uncertainty in this method is approximately a fac- 
tor of 2 (R. Maddalena 2002, private communica- 
tion). A value in the range 1 — 8 x 10^'^ cm~^ may 
therefore be considered typical for this cloud, with 
some preference for the lower end of that range. 

3.3.2. CO(l-O) Intensity 

From Figure 1 of Williams & Maddalena 
(1996), the total Iqo profile over the whole cloud 
integrates to w 7.1 K kms~^ The CO(l-O) map 
for the cloud integrated over the same velocity 
range of 16-38 km s~^ as the Hi is shown in Fig- 
ures 2 and 4 of Williams & Maddalena (1996), 
based on earlier data of Maddalena & Thaddeus 
(1985).^ From Figure 2 of Williams & Maddalena 
(1996) the CO(l-O) brightness ranges from 3-7 
contours with a few positions reaching 10 con- 
tours. A range of 3-8 contours or 6-16 K km s~^ 
therefore seems representative of the CO data, 
again with some preference for the lower end of 
that range. 

'Note that the I — b maps in Williams & Maddalena (1996) 
axe similar in shape but brighter by about a factor of 2 when 
compared to the earlier map of 7co in Fig- 2 of Maddalena 
& Thaddeus (1985). This is not likely to be an effect of 
differences in the velocity range over which the data have 
been integrated, since that range is smaller than the range 
of 15-40 km s~^used by Maddalena & Thaddeus. We use 
the more recent results of Williams & Maddalena. 



3.3.3. n and Go for G216 -2.5 

The range of values for Ico and A^(Hi) ob- 
tained above is plotted as a black box on our Fig- 
ure 2. In the context of our model, these values 
describe gas with a density range of 0.7 — 5 x 10^ 
cm~^ and an incident FUV flux of 0.1-7 of the 
standard value Go near the Sun, with likely val- 
ues being approximately 2 x 10^ cm~^ and 0.8, 
respectively. 

Although no independent measurements of the 
total volume density are available for G216 -2.5, 
the value we have obtained is in the range gen- 
erally indicated for GMCs in the Galaxy (n = 
n(Hl) + 2n(H2) « 2n(H2) - 2 x 50 cm'^, e.g., 
Blitz (1993)). As to the FUV flux, Williams & 
Maddalena (1996) estimated Go ~ 1 from what is 
known about the two nearby young stars identi- 
fied in their study as likely to be responsible for 
the photodissociation. Our model thus agrees well 
with the observations. 

3.4. Approximate Detection Limits for H i 
and CO (1-0) 

The practical detection limits for the Hi 21 cm 
line are w 5 and « 100 K km s~^ for single-dish 
survey and interferometer array imaging obser- 
vations, respectively, corresponding to 7V(Hl) « 
1 x 10^^ and 2 x 10^° cm^^. These are shown as 
the dashed (single-dish survey) and dotted (inter- 
ferometer array) lines for A^(Hl) on Figure 2. The 
corresponding values for the 2.6 mm CO line are 
« 1 K km s~^ (dashed line: single dish survey) 
and w 10 K km s^-'^ (dotted line: interferometer 
array). In all cases, this assumes that the gas 
clouds are resolved; otherwise these limits need 
to be increased further by the ratio of the beam 
area to the cloud area. 

Locating the detection limits on our Figure 2, 
we can conclude that the range in parameter space 
over which the combination of Hi and CO (1-0) 
can provide useful constraints on n and Go is given 
approximately by Go < 100, 10 cm~^ ^ n < 10^ 
cm~^. These boundaries are set by the decreasing 
sensitivity of the H I emission at high FUV flux lev- 
els (-/V(H i) increases only logarithmically for high 
Go at constant n) and by observational detection 
limits (/co disappears at low n owing to insuffi- 
cient excitation of the transition; H I disappears at 
high n as the gas stays mostly molecular). 
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3.4- 1- Working with Synthesis Imaging Data 

Many obKcrvcrs who report synthesis imaging 
(interferometer array) data in the hterature use 
units of Jy bearn"^ instead of K. The required 
conversion factors are (unfortunately) telescope 
dependent; Appendix C provides an approximate 
recipe for these conversions. 

3.5. Conclusions 

We have shown that a combination of observa- 
tions in the 21 cm line of Hi and the 2.6 mm line 
of CO can provide a useful diagnostic for physical 
conditions in distant PDRs. Both of these lines 
are readily observed in the Galaxy and in near- 
by galaxies using ground-based single-dish and 
interferometer-array radio telescopes. The useful 
range in parameter space for this combined radio 
line diagnostic is approximately bounded by Go 
< 100, 10 cm~^< n < 10^ cm~^, a range that 
covers a large fraction of the volume of the ISM in 
galaxies. The unique layered morphology of the 
emission in PDRs provides a means of identify- 
ing that part of the Hi and CO emission that is 
related when that morphology is observable, but 
confusion along the line of sight makes this sep- 
aration difficult in the Galaxy. This method will 
be especially useful for the interpretation of high- 
resolution Hi and CO synthesis imaging data in 
nearby galaxies where the FDR morphology can 
be more easily identified.^ Improvements in the 
model parameters and in the details of the com- 
putations are anticipated in the near future that 
will permit even more precise values of Go and n 
to be obtained. 

We are grateful to Bruce Draine and David Hol- 
lenbach for discussions on the physics of PDRs. 
David also provided very helpful comments on ear- 
lier drafts of this paper. The referee offered many 
useful suggestions that have also been incorporat- 
ed. 



However, in this "extragalactic" case, the effects of beam- 
smearing are severe and nmst be carefully considered in 
relating the observational parameters of the models to the 
measured brightnesses and profile velocity widths. 



A. COMPARISON OF ANALYTIC AND NUMERICAL MODELS FOR H i PRODUCTION 

IN PDRs 



The Hi column density in a PDR is calculated with the same physics used to determine the excitation 
of the H2 near-infrared fluorescence lines. A logarithmic form for the analytic solution to the equation for 
formation-destruction equilibrium was first given by Sternberg (1988) (sec also cq. A2 in Hollenbach et al. 
(1991)^ and eq. 3 in Goldschmidt & Sternberg (1995)). The model is a simple semi-infinite slab geometry 
in statistical equilibrium with FUV radiation incident on one side. The solution gives the steady state Hi 
column density along a line of sight perpendicular to the face of the slab as a function of x, the incident UV 
intensity scaling factor (see Appendix B), and the total volume density n of H nuclei. Sternberg's result is: 

(Al) 
where; 



D = the unattenuated H2 photodissociation 


R 


= the H2 formation rate coefficient 


rate in the average ISRF, 




on grain surfaces, 


a = the effective grain absorption cross section 


X 


= the incident UV intensity scaling 


per H nucleus in the FUV continuum, 




factor. 


A'' (Hi) = the Hi column density, 


n 


= the volume density of H nuclei. 



Equation Al has been developed using a simplified three-level model for the excitation of the H2 molecule 
and is applicable for low-density (n < 10^ cm~^), cold (T < 500 K), isothermal, and static conditions, and 
neglects contributions to A''(Hi) from ion chemistry and direct dissociation by cosmic rays. The quantity G 
here (not to be confused with Go used previously) is a dimensionless function of the effective grain absorption 
cross section a, the absorption self-shielding function /, and the column density of molecular hydrogen N2: 

G= CT/e-^-^^^dTV^. 
^0 

The function G becomes constant for large values of N2 due to self-shielding (Sternberg 1988). Using the 
parameter values in this equation adopted by Madden et al. (1993), we have: 

N{Ri) = 5 X 1020 X ln[l + (90x/n)] 

where n is in cm~^. This is a steady state model, with H2 continually forming from Hi on dust grain 
surfaces, and Hi continually forming from H2 by photodissociation. In order to compare this result with the 
"standard model" computations we need to relate x to the Go used by Kaufman et al. (1999); this is because 
Sternberg (1988) and Kaufman et al. (1999) use different normalisations for the FUV flux (see Appendix B). 
When distributed sources illuminate an FUV-opaquc PDR over 2tt sr, the conversion is x = Go/0.85 (see 
Footnote 7 in Hollenbach & Ticlcns (1999)), resulting in 9Gx/n = 106Go/n. 

We have fitted the analytic expression for iV(Hi) to the model computations in the range in which cosmic 
ray dissociation is not a major contributor, roughly for Gq > 1, n > 10 cm~^. The result is that no consistent 
improvement is obtained by using any value for the coefficient of Gq other than the value 106 deduced above, 
although a modest improvement is obtained by using a slightly larger value for the leading coefficient in 
the equation, 7.8 x lO^*^, corresponding to a value of 1.3 x 10"^^ cm^ for the effective grain absorption cross 
section. With these small adjustments, our final equation is: 

iV(Hi) = 7.8 X 10^0 x ln[l + (106Go/n)] cm'^ (A2) 

'David Hollenbach lias kindly pointed out to us that there is a typographical error in eq. Al of their paper; the right side should 
read Io{p/N^-^)n2e-^^ . 
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resulting in an value for the fit that is everywhere in excess of 99% over the fitted range of n and Gq. 

In Figure 3 we show values from equation A2 plotted as solid lines together with dotted contour lines 
from our standard numerical model. The agreement is generally good over much of the n-Go parameter 
space of interest here; differences occur mainly in the top left corner of the diagram, and at low values 
of FUV flux. In the top left corner the analytic formula under-predicts the amount of H I column density 
computed from the standard model by about 30% owing to H i production by ion chemistry reactions such 
as H++H2^ H+ + H, HCO++e- ^CO + H, and PAH-+H+ PAH + H (where PAH is a polycyclic 
aromatic hydrocarbon), which are important at high Gq and low n. At values of Gq < 1 the contours of 
A^(Hi) for the numerical model become vertical; this is because the standard model includes a low level of 
cosmic ray ionization (see main text) which contributes a small amount of Hi even for Go = 0. 

B. DEFINITIONS OF THE MEAN INTERSTELLAR RADIATION FIELD 

Kaufman et al. (1999) characterized the emission from photodissociation regions (PDRs) in terms of the 
"strength" of the FUV radiation field, Go, illuminating the medium at energies between 6 and 13.6 eV. 
Numerical values for Go were presented in units of the "Habing Field," following earlier PDR modeling by 
Tielens & Hollenbach (1985), who expressed Go in units of the "equivalent Habing (1968) flux of 1.6 x 10~^ 
ergs cm~^ s~^ appropriate to the average interstellar medium." However, other characterizations of the 
FUV background also appear in the literature for the "typical" and local ISM. Draine (1978) has developed 
an analytic fit to many of these, which he defined as the "standard UV background" from 5-13.6 eV. 

Although numerical values for the Habing (1968) and Drainc (1978) representations of the average inter- 
stellar background can be interconverted, Draine's Figure 3 makes it evident that they have different spectral 
shapes. Furthermore, the "strength" of the interstellar FUV field is described in the literature by several 
different quantities (e.g., energy density, photon intensity, energy flux) using a variety of units. Some authors 
use a single symbol to report normalized values of their field descriptor, while others specifically represent 
the normalized quantity by a ratio. Accordingly, the purpose of this appendix is to summarize the approach 
taken by Habing and Draine in developing their descriptions of the ISRF at ultraviolet wavelengths and to 
interrelate the representations of that field adopted by several authors whose work is germane to our paper. 

Habing (1968) calculated radiation energy density values, u\, at 1000, 1400, and 2200 A which were 
averaged "over a considerable area" in space, using bandwidths of 912 1040, 1300-1450, and 1900-2400 A, 
respectively. His recommended values for the averaged energy density {ux) were 

{ux) = 40 X 10"^^ ergs cm"^ A'^ at 1000 A (11.9 - 13.6 eV) 

= 50 X 10"^^ ergs cm"^ A'^ at 1400 A (8.6 - 9.5 eV) (Bl) 
= 30 X 10"^^ ergs cm"^ A'^ at 2200 A (5.2 - 6.5 eV) 

However, Habing's results depend strongly on assumptions made concerning the wavelength dependence of 
interstellar extinction and the presence of diffuse galactic light. Witt & Johnson (1973) noted that these are 
"two areas for which there existed little or no observational data at the time of his work." 

Draine (1978) converted these and other existing calculations of the interstellar UV energy density into a 
quantity that he defined as the angle-averaged photon flux 

F{E) = X^ux/iA-Kh^c) photons cm'^ s"^ sr"^ eV"^ (B2) 

where E is the energy, h represents Planck's constant, and c is the speed of light. The resulting values and 

a few observational data points arc plotted on his Figure 3, together with an analytical expression that is 
"in good agreement with all of the above results over the range 5-13.6 eV." That expression, 

Fd{E) = 1.658 X lifE - 2.152 x IQ^E'^ + 6.919 x IQ^E^ photons cm'^ s"^ sr"^ eV"^ (B3) 
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where E is measured in cV, was referred to by Drainc as the "standard UV background" in the ISM. Since 
it is computed on a "per steradian" basis, this quantity is actuaUy an angle-averaged photon intensity. 

We have used equation B2 to convert the Habing (1968) energy density values to angle- averaged pho- 
ton intensities (denoted Fh{E)) at the three wavelengths reported in his paper, and they agree well with 
the corresponding points on Figure 3 of Drainc (1978). We have also computed values for Draine's "stan- 
dard" UV background at the same wavelengths, using equation B3. At 1000 A, Fh{E) = 3.875 x 10^ 
photons cm-2 s'^ sr"! eV-\ Fd(^) = 6.627 x 10'"^ (in the same units), and Fh{E)/Fd{E) = 0.585. Thus, 
the photon intensity computed from Draine's fit at this wavelength is 1.710 times that found by converting 
Habing's calculated energy density to these units. This result is identical to the value computed by Draine & 
Bcrtoldi (1996) at 1000 A, and (to one decimal point) to the value given in footnote 7 of Hollcnbach & Tielens 
(1999). The same result could also be obtained by first converting Draine's photon intensity to an energy 
density for direct comparison to Habing's (1968) calculation. Draine & Bertoldi (1996) chose Habing's result 
at 1000 A as the basis for normalization because in "neutral regions, pumping of H2 is primarily effected by 
far-ultraviolet photons in the 1110-912 A range". The energy flux over 47r sr at 1000 A is F = A-kEFh{E) = 
1.2 X 10~^ ergs cm~^ s~^ = AttI, where / is the mean intensity. 

The values for these quantities at 1400 A are Fff(£') = 1.329 x 10^, FniE) = 2.611 x 10^, and 
Fh{E) / Fd{E) = 0.509. In this case, Draine's fit to the angle-averaged photon intensity is 1.965 times 
the converted Habing result, and the energy flux is F = AnEFHiE) = 2.1 x lO'^. At 2200 A, Fh{E) = 
3.094 X 10*^, Fd{E) = 3.748 x 10^, and Fh{E)/Fd{E) = 0.826. The "standard" Draine result is 1.211 times 
that calculated by Habing at this wavelength, and the energy flux F = 4ttEFh{E) = 2.0 x 10""^ . 

Draine's standard curve can be integrated from 6 to 13.6 eV to flnd the total photon intensity, i.e., 
1.547 X 10^ photons cm~^ sr~^. When this result is multiplied by the average energy per photon (i.e., 
1.407 X 10~^^ ergs) over that interval, the Draine flux over An sr is Fjj = 2.7 x 10"'^ ergs cm^^ s^^. If we 
assume that the photon intensity ratio found earlier in the 11.9-13.6 eV band, e.g., Fh{E) / Fd{E) = 0.585, 
applies over the entire FUV spectrum (i.e., 6-13.6 eV), then F^ = 1.6 x 10~^ ergs cm~^ s~^. This is the 
equivalent Habing flux. 

Since both ultraviolet normalization procedures are still used in the literature, confusion can arise when 
comparing descriptors of the ISRF. Kaufman et al. (1999) represent that field by an "incident FUV flux" , Go, 
that is already normalized by the equivalent Habing flux of the average ISM (see their Table 1). Normalization 
to the Habing field was also used by HoUenbach & Tielens (1999) and Wolfire et al. (1990). The use of 
this normalization in these models is purely historical, having been incorporated into early versions of the 
underlying code in support of PDR modeling begun in the 1980's (M. Wolfire 1999, private communication). 

In contrast, Sternberg (1988) and Sternberg & Dalgarno (1989, 1995) compute the "UV photon intensity" 
(i.e., x^v) by introducing a "UV intensity scaling factor" x which is defined to be unity for Draine's fit to 
the average interstellar field intensity 7^. In Sternberg (1988) and Sternberg & Dalgarno (1989), the Draine 
(1978) standard result is re-expressed as 

1.068 x 10-^ 1.719x10-2 6.853 x 10-1 , -2 1 tt 1 

h = ^2 ^ J3 photons cm ^ s ^ Hz \ (B4) 

where A is expressed in nanometers. When Sternberg's representation of Draine's standard background (i.e., 
eq. B4) is integrated from 6 to 13.6 eV, and the resulting value, 1.942 x 10* photons cm-^ s'^, is multiplied 
by the average energy per photon, the total FUV flux over An sr is Fs = 2.7 x 10^'^ ergs cm-^ s~^. Once 
again, we find that Fh ~ 1-6 x 10-^ provided the ratio Fh (E) / Fd (E) — 0.585 found earlier between 11.9 
and 13.6 eV (i.e., 912-1040 A) is assumed to apply over the entire FUV band. 

Equation B4 can be obtained from equation B3 by multiplying the latter by Aw sr and 4.136 x 10- eV 
Hz-i, provided E is replaced by hc/X = 1.240 x 10~^/\nm eV. Thus, 

= 1.6547r X 10-"fz5(WAnm). (B5) 
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We note that equation B4 represents a photon flux rather than the photon intensity, since it has been 
averaged over sohd angle. Sternberg & Dalgarno (1995) corrected this problem in their equation Al by 
multiplying equation B4 by l/47r, resulting in corrected units of photons cm~^ s~^ Hz~^ sr~^. When this 
correction is made, equation B5 becomes Ii, = 4.136 x 10~^^ Fu {he/ Xnm)- 

Madden et al. (1993) explicitly introduced field normalization through the ratio Xuv/Xo, where xuv 
is defined as the UV field intensity, although they also called it an FUV energy flux. Unfortunately, the 
subscript "UV" is dropped later in their paper where the ratio is shown as x/Xo, increasing the potential 
for equating the numerator (incorrectly) with Sternberg's parameter x upon casual inspection. The divisor, 
Xo = 2.0 X 10~* ergs cm~^ s~^ sr~^, was determined from the "UV field" results in the solar neighborhood 
presented by Draine (1978). Although not explicitly stated when xo was defined, their introduction identifies 
the band of interest as extending from 912 - 2000 A(i.e., 6.2 - 13.6 eV). Over 47r sr, the flux corresponding 
to the Madden et al. value of xo is Fm = 2.5 x 10~^ ergs cm~^ s~^, and Fjj = 1.5 x 10~^ if we continue to 
assume that Fh{E)/Fd{E) = 0.585. 

The key issue for this paper is how to relate the dimensionless fleld descriptor chosen by Kaufman et 
al. (1999), i.e., Gq, to the dimensionless quantity used by Sternberg (1988), X) or by Madden et al. (1993), 
Xuv/xo- Since x was normalized by Draine's standard background between 6 and 13.6 eV, and xo was 
computed by integrating the same curve over a similar energy range (the lower limits may differ by 0.2 eV), 
Xuv Ixo is equivalent to x- We illustrate this circumstance by comparing the flux computed earlier from 
Sternberg's representation of Draine's standard curve (i.e., from eq. B4), given by Fg = 2.7 x 10~^ ergs 
cm~^ s""'^, with the flux computed (using fewer signiflcant flgures) from the value of xo used by Madden et 
al. (1993). The latter is Fm = 2.5 x 10"*^ ergs cm~^ s~^. This conclusion is not affected by whether or not 
the Sternberg (1988) expression is multiplied by a solid angle factor (i.e., l/47r) provided his numerator is 
multiplied by the same factor. 

The dimensionless parameter Go = [47r / F{E)dE]/[4,TT J FH{E)dE] is a flux ratio that is comput- 
ed over the same spectral band as xO) where the denominator J FH{E)dE = 0.585 / Fn{E)dE in the 
highest energy band (i.e., 11.9 13.6 cV) addressed by Habing (1968). Accordingly, over that band Go 
= 1.710 / F{E)dE/ J F£){E)dE . Since xo (with units of ergs cm~^ s~^ sr~^) represents an intensity that is 
referenced to Draine's standard curve, it can be written as / Fjj{E)dE, whereby Go= 1.710x!7i//xo = l-710x. 
Inserting this into equation 12 of Madden ct al. (1993), the argument of the logarithm, i.e., [90 (xj/v/xo )/"■] + ! 
can be written as [53Go/n] + 1, in agreement with equation A2 when the latter is written to accommodate 
irradiation over 47r sr. 



C. RADIO BRIGHTNESS UNIT CONVERSIONS 

It is common practice to report results from imaging synthesis observations of radio spectral lines as 

contour diagrams of constant surface brightness using "J bcam"'^ km s"'^" as the unit.^ This unfortunate 
combination of units is telescope dependent through the use of the beam unit in the denominator. The 
conversion to telescope-independent units such as K km s~^ (which we use here in eqs. 1 and 2) requires 
knowledge of the beam area of the radio telescope in steradians, which is not well defined for any radio 
synthesis imaging telescope owing to the absence of information at short interferometer spacings. Fortunately 
(at least for now), most data processing systems perform an image restoration step that sets the final point 
spread function to be a two-dimensional Gaussian. This provides a definite (but not necessarily accurate) 
value for the two-dimensional integral over the beam, and thereby "solves" the unit conversion problem. If 
we take the FWHMs of this Gaussian to be 6i and 62, then the conversion for 21 cm H i image contours may 
be done with the relation 



^Or, quite incorrectly, as "J km s ^" (e.g. Sheth et al. (2002)). Another interesting (and equally erroneous) Vciriant is to give 
single dish spectra in units of "K beam"^" (e.g. Fig. 4 in Schinnerer & Scoville (2002)). 
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Tb(Hi) _ 6.057 X IQS 

dS/dn 61x62 ' ^ ' 

for Tb in K, dS/dQ in Jy beam~^, and 61 and ^2 in arcseconds. The corresponding relation for the 2.6 mm 
hne of CO(l-O) is: 

TbjCO) ^ 91.97 

with the same units as the previous equation. We emphasize again that these temperatures are well defined 
only for completely resolved sources. If the telescope beam does not resolve the source, the data can 
provide only a lower limit to the true brightness; such lower limits are usually called "main beam brightness 
temperatures" and labelled Tmb- 
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CO J=l-0 Intensity 



HI Column Density/10=' cm-^ 




Fig. 1. — a), (left panel) CO(l-O) emission at 2.6 mm from the surface of the model PDR as a function of 
the density n and incident FUV flux Go for the standard model parameters in Kaufman et al. (1999). The 
density n = ni +2n2 is the total density of H in atomic (ni) and molecular (712) form at the depth where the 
CO(l-O) line is formed. At this point the gas is essentially all H2. Contours of constant CO(l-O) emission 
Ico ai'c shown, labeled in both ergs cm~^ s"'^ sr"'^ and in CO-observer units of K kms~^. The PDR is 
assumed to have the dust-to-gas ratio 6 of the local ISM near the Sun. b). (right panel) Hi column density 
A''(Hl) in units of 10^^ cm~^ produced by photodissociation in the analytic approximation of eq. A2. This 
approximation is preferred for comparisons with observational data (see §2.2). 
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CO (J=l-0) Intensity and HI Column Density 




logio[n/cm-3] 

Fig. 2. Combined plot from the model computations of CO(l-O) emission Ico in Fig. la and A^(Hi) in 
Fig. lb. The contour labels, identical with Fig. 1, are omitted here for clarity. Shown are approximate 
observational limits for single dish (dashed lines) and interferometer imaging synthesis data {dotted lines). 
These limits are described in the text §3.4. The black "box" shows the approximate range allowed by the 
observations of G216 -2.5. 
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HI Column Density/lO^i cm-^ 
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Fig. 3. — Contours of constant H I column density 7V(H i) in units of 10^^ cm~^ in the PDR as a function of 
the density n and incident FUV flux Go for the standard model parameters in Kaufman et al. (1999) {dotted 
lines), and for the analytic approximation of eq. A2 {solid lines). The labeled contour values are for the 
numerical model; the contours for the analytic model are as in Fig. 1. 
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